System-bath entanglement in nanothermodynamics 
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We consider a quantum harmonic oscillator linearly coupled to a bath of harmonic oscillators 
and evaluate the degree of entanglement between system and bath using the negativity as an exact 
entanglement measure. We establish the existence of a critical temperature above which the system- 
bath negativity vanishes. Our results imply that system-bath entanglement is not responsible for 
the violation of the Clausius inequality observed in the low temperature/strong coupling regime 
[Phys. Rev. Lett. 85, 1799 (2000)], as the latter still occurs well above the critical temperature. 
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Thermodynamics is a cornerstone of modern physics. 
In the canonical description, it deals with systems that 
are in contact with a heat bath at constant temperature. 
An implicit assumption in the notion of thermal contact 
is that the system-bath interaction is vanishingly small, 
so that the total energy is simply the energy of the sys- 
tem plus that of the bath This approximation is 
well-justified in the limit of high temperatures, when the 
thermal energy of the system is much larger than the 
coupling energy. In the low-temperature limit, however, 
quantum-mechanical effects come into play and the in- 
teraction between the system and the bath cannot be 
neglected anymore Q . It has recently been shown 0, 0| 
that close to zero-temperature, a finite coupling strength 
leads to a violation of the familiar Clausius inequality, 
SQ < TdS, which states that the infinitesimal heat 5Q 
exchanged with the bath cannot exceed the product of 
the temperature T of the bath and the infinitesimal en- 
tropy change dS. Accordingly, work could be extracted 
from a single quantum heat bath. A scheme to exper- 
imentally verify this violation using nanoscale electrical 
circuits has been put forward in Ref. So far, the origin 
of the above deviation from standard thermodynamics in 
the quantum regime has been associated with entangle- 
ment between system and heat bath [H, 0, IS IS 0, Mam ■ 
However, this conjecture has never been tested explicitly. 

In this paper, we compute the system-bath entangle- 
ment using a microscopic model consisting of a quantum 
harmonic oscillator linearly coupled to a chain of har- 
monic oscillators (Rubin model) [T(| ■ Due to the linearity 
of the model, its dynamics is exactly solvable. Moreover, 
since its thermal states are Gaussian, we can rely on ex- 
act entanglement measures like the negativity An 
early study of entanglement between system and bath 
using the PPT criterion for 1 x N Gaussian modes [l2j 
has been performed in Ref. [HI]. In the following, we nu- 
merically investigate the dependence of the system-bath 
negativity on the coupling strength and the temperature, 
and compare it to the deviation from the Clausius in- 
equality. We find that the negativity vanishes above a 
critical temperature and that the Clausius inequality is 
violated even when system and bath are separable. Our 
results clearly demonstrate that the violation is not solely 



induced by system-reservoir entanglement. 

Thermodynamics in the quantum regime. The starting 
point of our analysis is the total Hamiltonian, 



H = H S + H E 



Hr 



(1) 



which describes a system S coupled to bath B via the 
interaction /. In classical thermodynamics, the inter- 
action energy is assumed to be small and therefore ne- 
glected. However, for very low temperatures, the ther- 
mal energy of the system is not necessarily much larger 
than the interaction energy and the coupling term can- 
not simply be discarded. This is particularly relevant 
for solid-state nanodeviccs, such as nanomcchanical os- 
cillators, which can be strongly coupled to their envi- 
ronment |14| . For finite coupling strength, the system 



Hamiltonian does not in general commute with the to- 
tal Hamiltonian, [Hs,H] ^ 0. Since there are then no 
joint eigenstates of H$ and H , a system strongly coupled 
to a zero-tcmperature bath is unavoidably in an excited 
state, in stark contrast to normal thermodynamics [l5l ]. 
Consequently, the system can perform positive work by 
relaxing to its ground state. It is important at this point 
to emphasize that the second law is not violated, since 
the produced work cannot exceed the work initially re- 
quired to couple system and bath 0, H, [l6| . 

In order to perform a quantitative study of the vio- 
lation of the Clausius inequality, we consider the Rubin 
model of a heavy harmonic oscillator coupled to a closed 
chain of harmonic oscillators [lj| • This model has played 
a major role in the microscopic understanding of Brow- 
nian motion [l7l ]. The total Hamiltonian of the Rubin 
model is of the form ([TJ) with: 



H, 



P 

2M 



(2) 



N r 



a=l 



2m + 2 



N-l 



J2 9(^+1 



H T = j - 



(q - xi) 2 + {x N - qf 



(3) 
(4) 
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For simplicity, the N oscillators of the chain have iden- 
tical mass m and frequency u>b- The respective mass 
and frequency of the system are denoted by M and u>s, 
and the coupling strength is given by the spring constant 



j/4. The Rubin model is fully characterized by 



the spectral density function (u> < ljr), 



J(lu) = — \Jlu 2 + UJ R UJ B \/ (bJR - u b )ur - ■ (5) 

In the limit of large system mass, M 3> m, and small 
bath frequency, ojb <C ojr, the Rubin model reduces 
to a Caldeira-Leggett-type model with Ohmic dissipa- 
tion, J{lo) = ^MloT 2 d /{lo 2 + T 2 D ), with friction coefficient 
7 = mu)ft/2M and Debye cutoff frequency Tjj = [13 • 
System and bath are furthermore supposed to be initially 
decoupled with a total density operator given by p(0) = 
Ps(0) ® p B and p B = exp(-(3H B )/Z B with (3 = (fcT)" 1 . 

The dynamics of the system is conveniently described 
using the phase-space representation. The Wigner func- 
tion of the system is found to satisfy the equation 0, [l8| , 

^W{ q , P ,t) = -^W + ^{[, P + M.h]W 



+ ~/MD pp (q,t) 



dp- 



■W 



dqdp 



(p„(q,t)w) . (6) 



In the high-temperature limit, the generalized diffusion 
coefficient D qp vanishes and Eq. ([6]) reduces to the usual 
Klein-Kramers equation [TtJ • The stationary solution of 
Eq. ([S]) is given by a quasi-Gibbs distribution, 



W(q,p) 



1 



2tt W T p T x 



cxp 



P 



Mu; 2 q 2 



2MT„ 



2T„ 



(7) 



with effective temperatures defined as T q = Mojg(q 2 ) and 
T p = (p 2 )/M. The stationary position and momentum 
quadratures can be calculated explicitly using Eq. ([5|) in 
the Ohmic regime and read 



T 



h \ - 



O 4=1 



(Ai-r D ) ip i + 



27T 
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E 
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- A,-) 



(8) 



(9) 



In the above equations, Ai are the characteristic fre- 
quencies of the damped harmonic oscillator and ij) de- 
notes the digamma function. In the limit of high tem- 
peratures/weak coupling, the two temperatures become 
equal, Tq = T p = T, and Eq. reduces to the usual 
Gibbs distribution. However, in the opposite limit of low 
temperatures/strong coupling, the system is squeezed by 
the coupling to the bath and T q 
ations from standard thermodynamics appear. 



-q < T p . As a result, devi- 
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FIG. 1: (Color online) Deviation from the Clausius equality, 
A = 8Q — TdS, for increasing values of the friction coefficient 
7: 7 = 0.3 (black narrow dashed), 7 = 0.6 (red solid) and 
7 = 0.9 (blue dashed), for the parameters M = 10, m = 1, 
us = 5, h = 1 and k = 1. 



We next focus on the violation of the Clausius inequal- 
ity in the low temperature/strong coupling regime. The 
internal energy U of the system is defined as the station- 
ary expectation value of its energy, that is, 



U = (H,) = 



(P 2 



(10) 



2M 2 

A quasistatic variation of a system parameter then yields 

dU = J J dqdpW(q,p) dH s {q,p) + 

J J dqd P H s (q,p)dW(q,p) . (11) 

The first term on the right hand side is identified with the 
work SW done on the system and the second term with 
the heat 5Q exchanged with the bath. In the following, 
we consider a variation of the mass M of the system. The 
change in heat is accordingly [3, [3] , 



5Q = 



(Mw 2 d(q 2 ) 



1 d(p 2 



V 2 dM 2M dM 



dAI. 



(12) 



On the other hand.the von Neumann entropy S of the 
system is given by [1, 3 , 



S = - tr{p s ln p s } 
= (v + -) In (v + 



(13) 



where ps is the stationary density oper ator of th e system 
corresponding to Eq. ([7]) and v = yj (q 2 )(p 2 ) jh. It is 
useful to introduce the quantity A = 8Q—TdS, where dS 
is defined as dS = OS/dMdM. The Clausius inequality 
then implies that A < 0, the equality being only achieved 
for quasistatic transformations [lj ■ Figure Q] shows the 
deviation A as a function of the temperature of the bath 
for different values of the friction coefficient 7. We clearly 
recognize that the Clausius inequality is violated, A > 0, 
for low temperatures, the violation being stronger the 
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larger the value of the coupling 7. On the other hand, 
with increasing temperature the deviation tends to zero, 
as expected for the quasistatic mass variation that we 
examine. It is commonly believed that this violation of 
the Clausius inequality in the quantum domain is due 
to the entanglement created between the system and the 
bath as a result of their interaction 0, 0, 0, E| 0, B 0] • 
We now show that this assumption is not correct. 

System-bath entanglement. Quantifying bipartite en- 
tanglement for general mixed states is a non-trivial task. 
However, for the special class of Gaussian states, such as 
the thermal oscillator states of the Rubin model, exact 
entanglement measures do exist. In the present investi- 
gation, we use the negativity to characterize system-bath 
entanglement. The negativity is defined as [ll|, 



P Ts H-l 



(14) 



where p is the total density operator of system plus bath 
and p Ts denotes the partial transposed with respect to 
the system. The trace norm of an operator A is given by 

The negativity is equal to the sum 

of the modulus of the negative eigenvalues of p Ts and 
therefore quantifies the degree of non-positivity of this 
operator. The negativity is zero for separable states and 
increases with increasing degree of entanglement. 

The Gaussian states of the Rubin model are fully char- 
acterized by the covariancc matrix, 
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FIG. 2: (Color online) System-bath negativity TV as a func- 
tion of the coupling constant 7 for different temperatures: 
T = 0.1 (red boxes), T = 1 (blue triangles) and T = 2 
(black diamonds), for the parameters M = 10, m = 1, wg = 5, 
ujb = 0.01, h — 1 and k = 1. 
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r ife w = -(oa+a^>a-(o) CT (a)c 



(15) 



where £ = (q,p, %i,Pi, ■ ■ ■ , xn,Pn) T is the vector of the 
combined positions and momenta of the system and bath 
oscillators. The symbol ( ) a represents the expectation 
value taken with respect to a density operator a. The 
trace norm \p Ts \ in Eq. (fl4|) can then be directly ex- 
pressed in terms of the symplectic eigenvalues i/j of the 
covariance matrix T(p Ts ) [llj : 



N+l 

n 



1 



for 
for 



> 1/2 
< 1/2 



(16) 



The symplectic eigenvalues of the matrix T(p Ts ) are the 
eigenvalues of the product QT(p Ts ), where is the sym- 
plectic matrix; they occur in conjugate pairs iuj, —ivj. 
For the case of (1 x TV) Gaussian modes that we con- 
sider, the negativity can be written in the simple form 

m, 



Af(p) 



\ 4f m in 




for v n 
for v n 



< 1/2 
> 1/2 



(17) 



with v mnl = min-jVj} the smallest symplectic eigenvalue. 

Numerical results. We have numerically computed the 
system-bath negativity, Eq. (|17p . and examined its de- 
pendence on the coupling constant 7 and bath temper- 
ature T. To this end, we have evaluated the full time 



FIG. 3: (Color online) System-bath negativity TV as a func- 
tion of the bath temperature T for different coupling con- 
stants: 7 = 0.3 (red boxes) 7 = 0.6 (black dots) and 7 = 0.9 
(blue triangles), for the same parameters as above. 



evolution of the coordinate vector £(t) and determined 
the covariance matrix (|15[) usin g the Williamson normal 
form of the Hamiltonian [2ll . |22| ; the system was taken 
to be initially in the ground state. The symplectic eigen- 
values Uj of the partially transposed covariance matrix 
were then computed in the stationary regime. The re- 
sults of our numerical calculation for a bath of N = 200 
oscillators are summarized in Figs. [2] and [3] 

Figure [5] shows the negativity as a function of the cou- 
pling constant 7 for various bath temperatures. As ex- 
pected, entanglement between system and bath increases 
with increasing coupling strength. However, for high 
temperatures, the existence of a critical coupling inten- 
sity 7 C below which the negativity is zero is clearly vis- 
ible. The presence of a threshold is further confirmed 
in Fig. [3] showing the temperature dependence of the 
negativity for different coupling constants. We observe 
that system-bath entanglement decreases with increas- 
ing temperature and completely vanishes above a critical 
temperature T c . The value of the critical temperature 
grows with growing coupling strength. 

It is instructive to compare our findings with the re- 
cent entanglement phase diagram proposed by Anders 
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FIG. 4: (Color online) Deviation from the Clausius equality 
A (red line) and system-bath negativity N (black dots) as 
a function of the bath temperature T for 7 = 0.6. All other 
parameters are those of Figs lH3l The violation of the Clausius 
equality is significant above the critical temperature at which 
the negativity abruptly drops to zero. 

and coworker for a closed chain of identical oscillators in 
thermal equilibrium [23|, [24| ; the Hamiltonian of the lat- 
ter is given by Eqs. © -(01) with M = m and lus = ^b- 
Entanglement between adjacent oscillators in the identi- 
cal chain was shown to vanish at a temperature kT c = 
or in our case, kT c = H , jM/2m. This esti- 
mate of the critical temperature, or cquivalently of the 
critical coupling, hr/ c = 2mkT/M, for a fixed tempera- 
ture, agrees very well with the values found for the Ru- 
bin model in Figs. M and H kT c ~ 5tvy for M = 10 and 
m = 1. We therefore obtain the important result that 
system-bath entanglement disappears at approximately 
the same temperature as entanglement within the bath. 
Violation of the Clausius inequality. The relationship 



between the violation of the Clausius inequality and the 
presence of system-bath entanglement can now be es- 
tablished by combining the results of Figs. Q] and [3J Fig- 
ure 3] unambiguously shows that the deviation, A > 0, 
from the Clausius inequality persists at temperatures way 
above the critical temperature at which the negativity 
vanishes. We can therefore conclude that the deviation 
from ordinary thermodynamics is not caused by the en- 
tanglement of the system with the bath, contrary to what 
has been claimed up to now. We note in addition that 
this discovery has also deep consequences for the theory 
of decoherence, where it is often believed that the loss 
of phase coherence originates from the entanglement of 
the system with the external reservoir (2f|. Again, this 
cannot be the case for temperatures above T c . We stress 
that the negativity can already vanish in the regime of 
low temperatures, kT c < huig, as illustrated in Fig. [3| 

To summarize, we have performed a detailed study of 
the entanglement between system and bath in the Rubin 
model of a heavy oscillator coupled to a chain of harmonic 
oscillators and showed that it vanishes above a critical 
temperature T c . We have further provided compelling ev- 
idence that the violation of the Clausius inequality which 
occurs in the regime of low temperature/strong coupling 
is not induced by system-bath entanglement, as it per- 
sists well above T c . 
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